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FOREWORD 


This  report  describes  an  effort  conducted  by  the  University 
of  Dayton  Department  of  Mathematics  as  a  subcontractor  to 
Universal  Technology  Corporation  under  Air  Force  contract  F33615- 
82-C-5089.  The  task  title  is  Statistical  Treatment  of 
Aramid/Epoxy  Stress-Rupture  Data  and  it  is  task  number  70.  The 
task  Mas  initiated  by  T.J.  Reinhart  through  AFWAL/Materi a1 s  Lab. 

The  work  reported  herein  was  performed  during  the  period  of 
15  Feb  1984  to  30  Apr  1984  by  the  author  Dr.  Gerald  Shaughnessy 
(University  of  Dayton).  The  report  was  released  by  the  author  in 
August  1984.  . 

The  author  would  like  to  thank  R.E.  Glaser  of  Lawrence 
Livermore  Laboratory  for  providing  the  computer  program  and 
output  which  made  the  calculations  in  this  report  possible. 
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1  Summary  of  stress-rupture  data  for  Aramid/Epoxy  strands  2 

2  Estimated  lifetime  of  Aramid/Epoxy  composite  under  sus-  7 

tained  tensile  loading.  (MLEL  calculation) 

3  Estimated  lifetime  of  Aramid/Epoxy  composite  under  sus-  8 
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4  Estimated  lifetime  of  Aramid/Epoxy  composite  under  sus-  9 

tained  tensile  loading.  (99%  confidence  interval) 


1.0  INTRODUCTION 


In  the  report  (4)  by  G.  Shaughnessy  and  J.  Snide,  the  experi¬ 
mental  stress-rupture  data  on  S-gl iss/apoxy  of  T.  T.  Chiao  (1)  was 
summarized  into  tabular. form  which  then  was  incorporated  into  the 
appropriate  AFSC  design  manuals.  This  report  will  do  a  simi^:r 
summary  of  the  stress-rupture  data  on  Aramid/epoxy  composite  :;trand 
data  collected  by  T.T.  Chiao  and  his  associates  at  the  Lawrence 
Livermore  National  Laboratories.  As  in  (4),  this  report  will  pro¬ 
vide  tables  of  estimates  of  quantile  levels  of  the  stress-rupture 
life  of  strands  subjected  to  sujtained  tensile  loading.  The  form 
of  the  data  that  was  collected  is  given  in  Table  1. 


TABLE  1 


SUMMARY  OF  STRESS-RUPTUkE  DATA 
FOR  ARAMID/EPOXY  STRANDS 


PERCENTAGE 

OF  FAILURE 
STRESS 

STRESS, KSI 

#STRANDS 

TESTED 

IEXACT 

FAILURE 

TIMES 

#GROUPED 

FAILURE 

TIMES 

fCENSORED 

TIMES 

90.0 

455.2 

101 

99 

2 

0 

87.0 

440.0 

100 

100 

0 

0 

84.0 

424.8 

103 

103 

0 

0 

80.0 

404.6 

100 

100 

0 

0 

70.0 

354.0 

49 

47 

2 

0 

60.0 

303.5 

50 

41 

9 

0 

50.0 

252.9 

50 

4 

0 

46 

2.0  ANALYSIS 


The  resulting  lifetimes  of  the  stress-rupture  tests  at  each  given 
stress  level  were  modeled  as  c  Weibull  distribution.  The  distribution 
function  for  a  Weibull  distribution  is 


F{t)  =  1  -  expC-(t/6)''] 


(1) 


where  F(t)  is  the  probability  that  the  lifetime  of  a  given  strand 
subject  to  a  stress  (o)  will  be  less  than  or  equal  to  t,  and  a  and B  are 
the  shape  and  scale  parameters  of  the  Weibull  model  to  be  estimated 
from  the  experimental  data.  These  parameters ,  a  and  8  of  the  Weibull 
model,  vary  with  the  given  stress  level  (o). 

The  basic  method  of  analysis  used  in  this  paper  to  estimate  the 
values  of  a  and  B of  the  Weibull  model  is  known  as  the  method  of  maximum 
likelihood.  In'  this  method,  the  joint  distribution  of  the  observed 
data  (called  the  likelihood  function)  is  regarded  as  a  function  of  the 

A  A  , 

parameters  Q and  B  .  The  values  a  and  B  that  maximized  the  likelihood 
function  are  referred  to  as  maximum  likelihood  estimates  of  the 
parameters  a  and  B  .  These  values  can  then  be  substituted  into  the 
reliability  equations  to  create  the  tables- 

As  stated  earlier,  the  parameters  o  and  B  vary  with  each  stress 
level.  The  following  two  relationships  were  used  to  describe  how  a 
and  B  varied  with  the  stress  level  for  the  data  being  analyzed. 

iTT)  •  9,  "  9,^  ^  9^0* 

*  * 

ln(8(o))  »  0  4.  00  4.00  +  00 


(2) 

(3) 


3 


9  are  constant  coefficents  that  are  to  be  determined.  Other 
n 

representations  of  the  relationship  between  j  cind  8  and  the  stress 
level  (a)  were . cons i dered  but  the  quadratic  '  odel  foroand  cubic  model 
for  8  fit  best  with  the  experimental  data  [2]. 

Using  these  representations  of  a  (q)  and  8  (a),  Eqs.  3  and  4  were 
then  substituted  into  the  likelihood  equation  which  was  then  maximized 
as  a  function  of  the  seven  variables  (9, 0,9,  0,9, 9, 6,). 

1  2  3  I*  S  S  T 

yv  A  A  A 

The  maximum  likelihood  estimates  (9  ,9  ,9  ,9  ,9  ,9  ,9  )  are 

.  l  2  3  I*  5  6  7 

then  obtained  by  setting  seven  partial  derivatives  of  the  likelihoocl, 
equations  equal  to  zero  and  solving  the  resulting  seven  equations  and 
seven  unknowns. 

Once  the  parameters  a  and  8  of  the  Weibull  model  have  been 
estimated  bv  the  above  procedure,  quantile  points  in  the  lifetime 
distribution  of  the  Aramid/Epoxy  composite  can  be  estimated.  A 
quantile  is  a  point  in  the  oistrlbutioh  that  has  a  given  percentage  of 
the  di stri buti Oil  less  than  it.  For  example,  the  0.1  quantile  in  a 
lifetime  distribution  is  a  time  such  that  the  probability  of  a  strand 
failing  before  that  time  is  0.1.  Using  the  maximum  likelihood 

A  ^ 

estimates  (a  and  8),  the  pth  quantile  in  the  lifetime  di  stri  Isuti  on  is 
obtained  by  setting 


(6) 


3.0  Results 


Table  2  was  generated  from  £q.  7  and  gives  the  maximum 
likelihood  estimated  life  (MLEL’  for  failure  probabilities  ranging  from 
0.0001  to  0.1  and  for  stress  levels  ranging  from  2G%  to  45%  of  the 
ultimate  flow  stress  (UTS).  ' 

For  example,  at  a  stress  level  of  30%  (UTS)  and  a  quantile 
value  of  0.001,  the  table  value  is  24.6  years.  This  means  that  the 
rraximum  likelihood  estimate  of  the  point  in  the  distribution  of 
lifetimes  under  which  there  is  a  probability  cf  failure  of  0.001  is 
24.6  years.  This  estimate  has  uncertainty  in  it  sirtce  it  is  baseo  on 
estimates  of  a  and  B  obtained  from  the  experimental  data.  Tables  3 
and  4  give  lower  confidence  bounds  at  95%  and  99%  levels  of 
confidence,  respectively,  for  quantiles  in  the  life  time  distribution 
under  which  the  probability  of  a  strand  failing  is  either  0.0001, 

0.001,  O.Cl,  or  0.1.  The  second  sat  of  tables  are  more  conservative 
and  reflect  the  size  of  possible  errors  in  the  estimation  procedures. 

At  the  95%  leve^  of  confidence,  the  values  in  the  table  were 
constructed  so  that  there  is  a  probability  of  0.95  that  the  true 
quantile  points  in  the  1  ifetime  distribution  are  not  lower  than  the 
table  values.  To  illustrate  further  the  use  of  the  tables,  suppose 
that  an  estiniata  of  the  0.001  quantile  point  in  the  lifetime  , 
distribution  of  a  strand  subject  at  a  stress  of  ,30%  (UTS)  is  desired. 
This  is  a  point  in  the  distribution  below  which  there  is  a  probability 
of  0.001  that  a  strand  will  fail.  Using  Table  2  for  the  MLEL,  the 
estimate  is  found  to  be  24.6  years.  Using  Table  3  for  the  lower  95% 
confidence  bound,  it  is  found  that  there  is  only  a  5%  chance  that  this 
0.001  quantile  is  below  3.5  years. 


4.0  CONCLUSIONS 


In  using  the  tables  in  thisreport  for  design  of  pressure  vessels 
for  Air  Force  app  1  i cat i on, i t  should  be  remembered  that  the  analysis  was 
based  solely  on  strand  data.  Also,  that  all  the  estimates  are  based  on 
extrapolation  from  the  experimental  data.  No  failure  data  was  observed 
below  50%  UTS.  Extrapolated  estimates  from  the  data  have  a  large  amoun‘ 
of  variability  in  them.  This  large  variability  resulted  in  lower  confi 
dence  bounds  fcr  quantiles  below  -C. 01-  that  are^  of  little  value  since  the 
variability  of  the  estimates  are  so  large. 

There  is  evidence  that  the  strand  data  is  conservative  for  design 
of  pressure  vessels  [3].  The  strands  tested  were  weakened  by  the  pre¬ 
sence  O'  f 1 uorescen t -  1 i gh t s  wh i ch  could  be  controlled  in  application. 
Also,  stress  rupture  data  on  vessels  have  indicated  that  strand  data 
is  conservative  in  redicting  vessel  strength  [3]. 


TABLt  2 

ESTIMATED  LIFETIME  OF  ARAM  1 0/ EPOXY  COMPOSITE 
JriDER  SUSTAINED  TENSILE  LOADING, 

MLEL  CALCULATIONS  (YEARS) 

(MAXIMUM  LIKELIHOOD  ESTIMATED  LIFE) 


QUANTILE  LEVEL 


5JFIBER 

STRESS 

STRESS 

(KSI) 

0.0001 
MLEL  (YRS) 

0.001 

MLEL  (YRS) 

0.01 

MLEL  (YRS) 

0.1 

MLEL( YRS) 

26 

131.5 

20.9 

44.5 

95.1 

205.8 

27 

136.6 

18.5 

37.8 

77.1 

159.4 

28 

141.6 

16.6 

32.4 

63.4 

125.5 

29 

146.7 

15.0 

28.1 

52.8 

100.4 

30 

151,7 

,  13.6 

24.6 

44.5 

81  .5 

31  ' 

156.8 . 

12.5 

21.8 

38.0 

_  '  67.1 

32 

161.9 

'11.5 

19.4 

32.3 

5  5.9 

33 

166.9 

10.6 

17.4 

28.5 

47.2 

34 

172.0 

9.8  . 

15.7 

25.0 

40.3 

35 

177.0 

9.1 

14.2 

22.2 

34,. 3 

36 

182.1 

8.5 

13.0 

19.8 

30.4 

37 

187.1 

7, .9 

11.9 

17.7 

26.7 

38 

192.2 

,7.4 

10-9 

16.0 

23.7 

39 

197.3 

6.9 

10.0 

14.5 

21.2 

40 

202.3 

6.5  , 

9.2 

13.2 

19.0 

4T 

207.4 

-  6,.0 

•8.5 

12.1  , 

17.2 

,  42 

212.4 

5.6 

7.3 

11. p 

,15.6 

43 

217.5 

5.1 

7.2 

10.1 

14.3 

44 

222.6 

4.7 

6.6 

9.3 

13.1 

45 

227.6 

4.3 

6.1 

8.5 

12.0 

TABLE  3 

ESTIMATED  LIFETIME  ,0F  ARAMID/EPOXV  COMPOSITE 
UNDE‘5  SUSTAINED  TENSILE  LOADING 

LOWER  BOUND  CALCULATIONS  (YEARS) 

95%  CONFIDENCE  INTERVAL 


QUANTILE 

LEVEL 

%FIBER 

STRESS 

STRESS 

(KSI) 

0.0001 

LB  (YRS) 

0.001 
-LB  (YRS) 

0.01 

LB  (YRS) 

0.1 

LB  (YRS) 

26 

131.5 

1.9 

3.5 

10.4 

24.9 

27 

136.6 

1.9 

3.5 

9.9 

22.6 

28 

141.6 

.  1.9 

3.5 

9.4 

20.7 

29 

146.7 

1.9 

3.5 

8.9 

19.1 

30 

151.7 

1.9 

3.5 

8.6 

17.7 

31 

156.8 

1-9 

•3.5 

8.3 

16.6 

32 

161.9 

1.9 

3.5 

8.0 

15.6 

33 

166.9 

I  .9 

3.5 

7.8  ■ 

14.7 

34 

172.0 

1.9 

3.5 

7.6 

14.0 

35 

177.0 

1.9 

3.5 

■  7.4 

13.4 

36 

182.1 

1.9 

3.5 

7.3 

12.8 

37 

187.1 

1.9 

3.5 

7.1 

12.3 

38 

1 9,2 . 2 

1.9 

3.5 

7.0 

1 1.8 

39 

19  7.3 

1.9 

3.5' 

6.3 

11.3 

40 

202.3 

1.9 

3.5 

6.7 

10.9 

41 

207.4 

1.9 

3.5 

6.5 

10.5  . 

42 

212.4 

1.9 

3.5 

6.3 

10.1 

43 

217.5 

1.9 

3.5 

6*1 

9.7 

44 

222.6 

1.9  , 

3.5 

5.9 

9.3 

45 

227.6 

1.9 

3.5 

5.7 

8.8 

8 


TABLE  4 


ESTIMATED  LIFETIME  OF  ARAMID/EPOXY  COMPOSITE 
'JNDER  SUSTAINED  TENSILE  LOADING 

LOWER  BOUND  CALCULATIONS  (YEARS) 

99%  CONFIDENCE  INTERVAL 


QUANTILE 

LEVEL 

*FIBER 

STRESS 

STRESS 

(KSI) 

0.0001 

LB  (YRS) 

o.oo; 

LB  (YRS) 

O.'Ol 

LB  (VRS) 

0.1 

LB  (YRS 

26 

131.5 

1.4 

2.6 

5.0 

10.4 

27 

136.6 

1.4 

2.6 

5.0 

10.1 

28 

141.6 

1.4 

2.6 

5.0 

9.8 

29 

146.7 

1.4 

2.6 

5.0 

9.6 

30 

151.7 

1.4 

2.6 

5.0 

9.4 

31 

156.3 

1.4 

2.6 

5 . 0 

9.3 

32 

161.9 

1.4 

2.6 

5.0 

9.2 

33 

166.9 

1.4 

2.6 

’  5.0 

9.1 

34  • 

172.0 

1.4 

2.6 

5.0 

9.0 

35 

177.0 

1  .4 

2.6 

5.0 

9.0 

36 

182.1 

■1.4 

2.6 

5.0 

8.9 

37 

187.1 

1.4 

2.6 

5.0 

8.9 

38 

192.2 

1.4 

2.6 

5.0 

8.8 

39 

197.3 

1.4 

2.6 

5.0 

8.7 

40 

202.3 

1.4 

2.6 

,5.0 

8.7 

41 

207.4 

1.4 

2.6 

5.0 

8.6 

42 

212.4 

1.4 

2.6 

5.0 

8.4 

43 

217. S 

1.4 

2.6 

5.0 

8.3 

44 

222.6 

1.4 

2.6 

4.9 

8.1 
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LIST  OF  SYMBOLS.  ABBREVIATIONS,  ACRONYMS 


MLEL  Maximum  likelihood  estimated  life. 

F(t)  Probability  that  the  lifetime  of  a  given  strand  subject 

to  a  stress  will  be  less  than  or  equal  to  the  time  the 
stress  is  appl ied. 

t  Time . 

alpha  Shape  parameter  from  Weibull  model, 

beta  Scale  parameter  from  Weibull  model, 

sigma  Stress, 

theta! -theta?  Constants, 

p  A  given  probability  F(t). 

UTS  Ultimate  flow  stress, 

ksi  Kips  per  square  Inch. 
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